The B and Bg Meson Decay Constants from Lattice QCD 
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We present a new determination of the B and Ba meson decay constants using NRQCD b-quarks, 
HISQ light and strange valence quarks and the MILC collaboration Nf = 2 + 1 lattices. The 
new calculations improve on HPQCD's earlier work with NRQCD 6-quarks by replacing AsqTad 
with HISQ valence quarks, by including a more chiral MILC fine ensemble in the analysis, and 
by employing better tuned quark masses and overall scale. We find /s ~ 0.191(9)GeV, fs, = 
0.227(10)GeV and fBs/fs = 1.188(18). Combining the new value for fBs/fs with a recent very 
precise determination of the Bs meson decay constant based on HISQ 6-quarks, fs^ = 0.225(4)GeV, 
leads to /b = 0.189(4)GeV. With errors of just 2.1% this represents the most precise /s available 
today. 
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I. INTRODUCTION 

Precision electroweak data gathered at the B facto- 
ries, the Tevatron and at LHCb is allowing particle physi- 
cists to carry out stringent tests of the Standard Model 
(SM) and search for hints of New Physics (NP). Sev- 
eral groups, for instance, are studying global fits to the 
Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle 
(UT) and checking whether various combinations of con- 
straints coming from experiment and theory can be ac- 
commodated consistently with each other [l-y] . In recent 
years some tensions at the 2-3 a level within the SM have 
emerged from these studies and it will be very interesting 
to see whether future improvements in experimental and 
theory inputs will remove these tensions or conversely 
elevate them to serious hints of NP. 

Lattice QCD is playing an important role in UT anal- 
yses, providing crucial inputs such as ex, Bb^, £, = 
Jb, Bb^I fs^Bs, Jb and information on semileptonic 
form factors Q. To make progress in resolving the ten- 
sions in UT analyses it is imperative to reduce the errors 
in current lattice results. In reference [l| the B meson 
decay constant /b is not used as an input for the global 
fits but becomes instead one of the fit outputs /]/**^ 
This is then compared with the SM (i.e. Lattice 

QCD) value f'^^^^ to check for consistency. The au- 
thors of reference [l| experiment with dropping different 
processes in their global fits and study how this affects 
/b^*'' and when Jb^*^ agrees best with fl^^^\ Using 
this fit-comparison procedure, the authors attempt to de- 
termine the dominant source of deviations from the SM, 



e.g. whether it is coming from B ipKs (sin2/3), Bs 
and Bd mixing, Kaon mixing [ck) or B — > tv. Needless 



to say this 



j(QCD) QQjjj^pj^j-jgQj^ method requires 



knowing f^'^^^ as accurately as possible. In this article 
we significantly reduce errors in 

fiQCD)_ ^^^^ reduced 
errors, the B meson decay constant will hopefully help 
further constrain UT analyses in the future. 

In the next section we introduce the lattice setup and 
explain how the bottom and strange quark masses were 
fixed in our lattice actions. Section III discusses operator 
matching between heavy-light currents in full continuum 
QCD and in the lattice theory. We describe two-point 
correlators and the smearings employed. In section IV 
we present our fitting strategies to the two-point corre- 
lators and describe how the extracted amplitudes lead 
to the hadronic matrix elements relevant for determining 
decay constants. This section also includes summary ta- 

/s, ^/MsT, * = JbVM^, and 



bles of fit results for $s 
their ratios for each of the 6 MILC ensembles that we 
work with. Then in section V we explain how simulta- 
neous extrapolations to the continuum and chiral limits 
(physical point) were carried out. Section VI discusses 
results at the physical point and the error budget and we 
conclude with a summary in section VII. For the rest of 

this article we omit the " QCD" 



■ AQCD) 



II. THE LATTICE SETUP AND TUNING OF 
BARE QUARK MASSES 

HPQCD's previous work on B and Bg meson de- 
cay constants with NRQCD 6-quarks used AsqTad light 



and strange quarks |5| 
Nf = 2 + 1 lattices 
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It utilized the MILC AsqTad 
In the present work we re- 
place the AsqTad valence quarks by their Highly Im- 
proved Staggered Quark (HISQ) [tJ counterparts thereby 
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TABLE I: Simulation details on three "coarse" and three 
"fine" MILC ensembles. 



Set 


ri/a 


mi/nia (sea) 


Nconf 


Ntsrc 


X Nt 


CI 


2.647 


0.005/0.050 


1200 


2 


24^ X 64 


C2 


2.618 


0.010/0.050 


1200 


2 


20^ X 64 


C3 


2.644 


0.020/0.050 


600 


2 


20^ X 64 


FO 


3.695 


0.0031/0.031 


600 


4 


40^ X 96 


Fl 


3.699 


0.0062/0.031 


1200 


4 


28^ X 96 


F2 


3.712 


0.0124/0.031 


600 


4 


28^ X 96 



reducing the dominant discretization errors coming from 
staggered taste breaking by roughly a factor of three. 
Details of the MILC ensembles employed here are given 
in Table L There is considerable overlap between the 
MILC ensembles used in the present article and in Q . In 
an additional coarse ensemble with sea quark masses 
rrii/nis — 0.007/0.05 was employed. Here we have added 
instead a third, more-chiral fine ensemble, the 40*^ x 96 
Set FO with mi/ms = 0.0031/0.031. 

For the 6-quarks in our simulations we use the same 
NRQCD action employed in ^] . Since the publication of 
the HPQCD collaboration has updated the value of 
the scale parameter ri to ri = 0.3133(23)fm [s*], and this 
necessitated a retuning of all quark masses including the 
bare 6-quark mass aMf, for all MILC ensembles in Table 
1. To fix aMt, we use the spin averaged T mass. One 
calculates. 



with 



2AE„ 



AEp^E{p)-E{0), (2) 



and compares with the experimental value (adjusted 
for the absence of electromagnetic, annihilation and sea 
charm quark effects in our simulations) of 9.450(4)GeV 
0. Results from this tuning are shown in Fig.l. Errors 
in the data points include statistical and ri/a errors. One 
sees that these are much smaller than the 0.7% error in 
the absolute physical value of ri . To achieve small statis- 
tical errors in Mkm it was crucial to employ random wall 
sources for the NRQCD 6-quark propagators. In Fig. 2 we 
show a comparison between local and random wall results 
for T correlators with momentum 27r/(aL). Most of the 
tuning of aM^ was carried out with momentum 2t: /{aL) 
for ensembles CI, C2, C3, Fl and F2, and with momen- 
tum 47r/(aL) for ensemble FO. However, we have checked 
on one ensemble that consistent Mkin values result from 
higher (but not too large) momenta as well. For instance 
on C2 with aMt, = 2.8 (slightly larger than the actual 
physical 6-quark mass) one finds aMkin{^ Si) = 5.933(15) 
for momentum 27r/(aL) and aMkini'^ Si) — 5.941(15) for 
momentum 47r/(ai). 
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FIG. 1: Tuning of the &-quark mass via the spin averaged 
T mass. 9.450GeV corresponds to the experimental value ad- 
justed for lack of electromagnetic, annihilation and sea charm 
quark effects in the simulations. 



The s-quark mass was tuned to the (fictitious) r]s mass 
of 0.6858(40)GeV 0. Fig.3 shows results for this tun- 
ing. All but the Set FO point (most chiral point on plot) 
were fixed already in [10] . Having fixed the bottom and 
strange quark masses on each ensemble one can investi- 
(1) gate the mass combination AIb, — The leading 

dependence on the heavy quark mass cancels in this dif- 
ference, so one is testing how well the lattice actions are 
simulating QCD boundstate dynamics. Results for this 
mass difference are shown in Fig. 4. Within the ri scale 
error and additional ^lOMeV uncertainty from relativis- 
tic corrections to M^j^ one sees agreement with experi- 



ment in the continuum limit. 

Table II summarizes the valence quark masses used in 
this article. We include the HISQ valence charm quark 
masses for each ensemble, since these provide a conve- 
nient scale in the chiral/continuum extrapolations of sec- 
tion V. The charm quark masses were fixed by tuning to 
the 77c mass. The light HISQ valence quark mass mi 
is chosen so that mi{valence)/ms(valence) is close to 

mi{sea)/ml'^XtqTad^ where m^^'j'^^^^^ corresponds to the 
physical AsqTad strange quark mass. As a final consis- 
tency check of our lattice setup, we have looked at the 
Bg — B mass difference. This is shown in Fig. 5 . 



III. OPERATOR MATCHING AND RELEVANT 
CORRELATORS 



Decay constants Jb^ are determined by calculating the 
matrix element of the heavy-light axial vector current 
between the Bq meson and hadronic vacuum states. For 
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FIG. 2: Comparison of relative errors in the finite momentum FIG. 4: The mass difi^erence Mb — M -/2. 

T correlator with and without random wall sources. This 
example is taken from the C2 ensemble. 
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FIG. 3; Tuning of the strange quark mass via the fictitious 
rjs meson. 



TABLE II: Valence quark masses 



Set 


ami 


arris 


anic 


aMt 


CI 


0.0070 


0.0489 


0.6207 


2.650 


C2 


0.0123 


0.0492 


0.6300 


2.688 


C3 


0.0246 


0.0491 


0.6235 


2.650 


FO 


0.00339 


0.0339 


0.4130 


1.832 


Fl 


0.00674 


0.0337 


0.4130 


1.832 


F2 


0.0135 


0.0336 


0.4120 


1.826 



FIG. 5: The Bs - B mass difference AM versus the light 
valence quark mass. 



the temporal component in the Bq restframe one has, 

{0\Ao\Bq)QCD=MB,fB,. (3) 

Simulations are carried out with effective lattice theory 
currents, 



-1 



(4) 
(5) 



4'\^) = ^^^^-^^oToVf-Q, (6) 
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where ^'g is the HISQ action Ught or strange quark field 
(in its four component "naive fermion" form) and the 
heavy quark field with the upper two components given 
by the two-component NRQCD fields and the lower two 
components set equal to zero. We have matched these 
effective theory currents to Aq in full QCD at one-loop 
through order as, ^'^'^^ , aas, as ^'^m" ■ Details 

of the matching of NRQCD/HISQ currents will be pre- 
sented in a separate publication [ll|. The calculations 
follow the strategy developed in |l2| and employed for 
NRQCD/AsqTad currents in [H]. One finds, 

{A^)qcd = (l + a,po)(Jr) + (l + ".Pi)(^o''''"''> 
+a.p2(jf '^"'), (7) 

Here poi Pii P2 and Cio are the one-loop matching coeffi- 
cients. 

We use smeared heavy-light bilinears to represent the 
Bq mesons. For instance, we create a meson at time to 
via, 

$„(x,to) = ^'^Q{xi,tQ)aa{xi ~x)Tsc^q{x,to), (9) 

with = 75- For the smearing functions (Ja{xi — x) 
we use a 5-function local smearing (a = 1) or Gaussian 
smearings oc /(2''o)for two different widths ro and 

normalized to one {a — 2,3). We then calculate a 3 x 
3 matrix of zero momentum meson correlators with all 
combinations of source and sink smearings, 

X y 

with V — L^. We use Gaussian widths in lattice units 
of size To = 3 or 5 on coarse ensembles and tq = 4 or 7 
on the fine ensembles. In addition to this matrix of B 
correlators we also need correlators with at the source 
and Jq at the sink for i — 0,1, 2, 

c'aitM ^^J2J2(Jo\y,t)^.is,to))- (11) 

X y 

Since 70 ^'q — it turns out that, 

(12) 

Furthermore for zero momentum correlators one can 
show that, 

Cf^Ci\ (13) 

so only the three C^^, a = 1, 2, 3, are required in addition 
to the 3x3 matrix C^^. 

The spatial sums ^ - in (ITUl) and pTt are done at the 
sink, and so can be handled very easily. We implement 



the J2s sums at the source via random wall sources. This 
is described for instance in reference Q. Here we give 
some of the explicit formulas. In terms of quark propa- 
gators for the and fields ea. pU)) becomes (we set 
to = for simplicity), 

X y Xl yi 

{tr{GQ{yi - xi,t)aa{xi - x)TscGq{x- y,-t) 

^skcrpim ~ y)}) 
-^EEEE(^-{Gq(27i--1'0 

X y Xl yi 

aa{xi - x)TsclbG\{y- x,t)^^T skfypiyi - y)})- 

(14) 

We set, 

GQ"^°'\yi.,x,t) = ^GQ{yi - xi,t)aaixi - x), (15) 

Xl 

and recall the relation between the naive HISQ prop- 
agator Gq{y — x) and the one component HISQ quark 
propagator G^{y — x) [l3 |. 

Gqiy-x)^n{y)G^{y-x)n^x), (16) 

or equivalently, 

Gl{y-x)^n{x)[n{y)G^{y-x)]\ (17) 

with, 

17(x)^7o°7r72^73^ (18) 
Setting Fsc = Tsk = 75 one has, 

c^At) = 

^EEE(^-{[gS'""^(j^1'-'*m-)] 

X y yi 

[n{y)G^{y~x,t)]^af,{y,-y)}). (19) 

We introduce a random U(l) field £,{x) at each spatial site 
of the source time slice (in practice we employ separate 
U(l) fields for each color but suppress this index in the 
formulas given below) and replace, 

^E--^EE^(-^)^^(-^'). (20) 

X X x' 

Ea. (fTO|) becomes, 

Clc^t) = 

^llllT.T.^t'{ [Gg""^ (yi' tMx)ai)] 

X x' y yi 

p{y)G^{y-x',mx')]^ap{m-y)]). 

(21) 
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An even more concise expression can be obtained if one 
defines, 

X 



(22) 



and, 



^^Y.^iy)G^{y-x',t)ax'). 



(23) 



This leads to, 

y yi 

{tr { [G^,^™"'"") (yi , t)] [G(™) {y, t)] V, (yi - y) }) . 

(24) 

Equations (|22l) and (1^51) tell us that we should create 
NRQCD propagators with source, 

5Gg(fi) = -^E'^"(^i-^M^)^(^)' (25) 



and HISQ propagators with source, 

sc,{x') = -^a^') 



(26) 



The double sum in ([M)) is carried out via Fast Fourier 
Transforms. 



IV. FITS AND DATA ANALYSIS 

The 3x3 matrix of correlators C^^ of eq. (|10l) and the 
Ci' of eq.dni) for i = 1 can be combined into a 4 x 3 
matrix of correlators CjB^a with CjB^a = for a,/3 = 

1,2,3 and Cp=4^a=i,2,3 = C';^!;^ 2,3- Various subsets of 
these correlators are then fit simultaneously to the form. 



N-1 



N-1 



3=0 



k=0 



(27) 



to extract the ground state energy Eq and amplitudes 
&Q. The hadronic matrix elements appearing in ([7]) are 
related to the amplitudes &q as. 



a' (0|J^°^|B„ 



2A/B,a &S 



0=1 



(28) 



1x1 (loc-loc) 1 


■ 1 


2x2 (loc-sml) 




2x2 (loc-sm2) 


i-e-i 


3x3 (loc-sml-sm2) 




3x2 (loc-sml-jl) 


-e- 


3x2 (loc-sm2-jl) 




simultaneous and B 





0.53 



0.535 



0.54 
aE(Bs) 



0.545 



0.55 



FIG. 6: Examples of results from different matrix fits for the 
Bs meson energy in lattice units. These fit results are taken 
from the C2 ensemble and used = iV = 7 in ea.(|27)l. 



and 



(0|J^ 



4^\b,) = J2MBah 



(29) 



The factors of yJ^Ms^a come about due to differences 
in normalization of states in the effective lattice theory 
compared to the standard relativistic normalization of 
states. 

We have investigated fits to various subsets of corre- 
lators (submatrices) taken from the full 4x3 matrix of 
12 correlators. For each correlator we fit data between 



t ty] 



and t 



on coarse ensembles and tr 



4 and 
8 and tm 



16 

24 



on the fine ensembles. In Fig. 6 we show results for the 
Bs energy in lattice units, qEb^, from fits to ensemble 
C2. One sees a large improvement upon going from a 
fit to a single local-local (a,/? — 1) correlator to a 2 x 2 
matrix of correlators {a, (3 — 1,2 or a, (3 — 1,3). There 
appears to be little further improvement when one goes 
to 3 X 3 matrices. Our final fit results are taken from 
3x2 matrix fits with a = 1,3 and f3 = 1,3,4. We do 
simultaneously a 3 x 2 fit to _B correlators together with 
a 3 X 2 fit to i?s correlators. This allows us to get ratios 
such as Jb, \/Mb^/ JbVMb and mass differences such as 
AI B^ — Mb in a single fit with correctly correlated errors, 
in addition to the separate quantities /b and Jb^ ■ 

In all our fits we use Bayesian methods [l5| and work 
with fixed tmin and t^ax while increasing the number 



of exponentials N and N in eq. l(27|) until fit results in- 
cluding errors and chisquares/dof have saturated. Fig. 7 
shows fit results for the B meson amplitde 6j on ensem- 
ble CI versus N (which we also set equal to N). One 
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TABLE IV: Bs and B mass difference AM, and ratios 
/>I>'-''^ and ^s/O, with <1> = /s, y^Ms^. Errors include 



-I- 



statistical and fitting errors. 



Set 


AM [GeV] 






CI 


0.0648(22) 


1.1311(90) 


1.1324(89) 


C2 


0.0577(18) 


1.1132(65) 


1.1143(64) 


C3 


0.0413(20) 


1.0772(81) 


1.0775(80) 


FO 


0.0717(29) 


1.1508(123) 


1.1516(121) 


Fl 


0.0614(20) 


1.1223(75) 


1.1234(73) 


F2 


0.0478(13) 


1.0889(51) 


1.0896(50) 



CHIRAL AND CONTINUUM 
EXTRAPOLATION 



FIG. 7: Fit results for the B meson amplitude fej on ensem- 
ble CI versus the number of exponentials N = N = Nexp- 
Simultaneous 3x2 matrix fits were carried out to both B 
and Bs meson correlators at the same time. 



TABLE III: $ = JbVM^ and $a = /s, a/A/s^ in lattice 

units. The lowest order results and are also shown. 
Errors include statistical and fitting errors. 



Set 


a3/2${0) 








CI 


0.2394(18) 


0.2214(16) 


0.2708(13) 


0.2508(12) 


C2 


0.2498(18) 


0.2313(17) 


0.2780(9) 


0.2577(8) 


C3 


0.2545(28) 


0.2356(26) 


0.2742(14) 


0.2539(13) 


FO 


0.1431(16) 


0.1293(14) 


0.1647(8) 


0.1489(7) 


Fl 


0.1483(12) 


0.1340(11) 


0.1664(6) 


0.1506(5) 


F2 


0.1520(12) 


0.1375(10) 


0.1656(7) 


0.1498(6) 



sees that things have stabilized by = 4. In Table 

,3/2 , 



III we collect fit res ults f or a^/^$ = a^^'''fB\/MB and 
o?'^'^ fs,^ Mb^ for the six ensembles. The 



quantities ^'^"^ and are analogous results if only the 
1 X ( Jo^"*) contribution is included on the RHS of eq.©, 
i.e. if one drops all one-loop and 1/M current corrections. 
In Table IV we summarize results for the mass differ- 
ence AM = Mb, - Mb in GeV's and the ratios <^s/^ 
and <^^s'' /<^^^\ Fig.5 illustrates the results for AM. For 
the ratios one sees good agreement between ^s/^ and 
<I>^°^/(I>(o) indicating complete lack of sensitivity to 0{a) 
or 0{1/M) current corrections in this ratio. 



In this section we describe how we extrapolate $, $s 
and ^s/^ given in Tables III and IV to the physical 
limit, a 0, mi/nis {mi /nis) physical = 1/27.4 and 
ms/mc -> (ms/mc) physical = 1/11.85. We fit $ and $s 
to the general form, 

= $0 (1 + Sfq + [analytic]) (1 -f- [discret.]), (30) 

where Sfq includes the chiral logarithm terms. Explicit 
expressions, taken from the literature [l6l, [l7j. are given 
in the appendix. Most of our extrapolations employed 
formulas for 6fq at one-loop order in chiral perturbation 
theory (ChPT) and at lowest order in 1/M. We have 
also included some 1/M corrections such as effects of 
the B* - Bq hyperfine splitting as discussed in [l^l- For 
the [analytic] terms we use powers of m^aience/'mc and 
msea/mc, where nric is the bare charm quark mass (see 
Table II) fixed for each ensemble through the 7]c mass, nric 
is a convenient scale to use since ratios such as ung/mc 
or mi/mc are equal to the corresponding ratio of MS 
masses and are furthermore scale independent (up to dis- 
cretization corrections). For the [discret] terms in ((30)) 
we employ powers of (a/ri)^. We allow for the expan- 
sion coefficients to be themselves functions of aMb and/or 
arUq to take into account that we are dealing with an ef- 
fective NRQCD theory for the 6-quarks and with taste 
breaking splittings in staggered meson masses. In prac- 
tice we find that adding such effects into [discret] has 
minimal effect on the final value for decay constants in 
the physical limit. 

Fig. 8 shows chiral/continuum extrapolation results for 
$ b ■ We show extrapolations for our basic ansatz with, 

[analytic] — 

j3o{2mu + ms)/mc + Pimq/mc + Mmq/rric)^, 

(31) 

where mu{mq) denotes the sea(valence) light quark mass, 
[discret.] = co{a/ri)'^ + ci(a/ri)^, (32) 
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FIG. 8: Chiral/continuum extrapolation of $ = fs^/ Mb- 
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FIG. 9: Tests of /s at the physical point. The left most 
magenta point is the "basic ansatz" result. The remaining 
points refer to results when the basic ansatz was modified in 
several ways as explained in the text. 



and using (A. 7) from the appendix for Sfq. The /dof 
for this fit was 0.24. We have checked the stability of 
our extrapolations by modifying the basic ansatz in the 
following way: 



1. dropping the (32 term in (j31[) : 



FIG. 11: Chiral/continuum extrapolation of ^s/^. 

2. adding more {mq/mcY^ terms with n > 2; 

3. dropping the ci term in (|32p : 

4. adding more (a/ri)" terms with n > 4; 

5. making the coefficients depend on ruq i.e. 
power series in mq/mc] 
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TABLE V: Error budget 



Source 




fB 


fsJfB 




(%) 


(%) 


(%) 


statistical 


0.6 


1.2 


1.0 


scale r^^^ 


1.1 


1.1 




continuum extrap. 


0.9 


0.9 


0.9 


chiral extrap. & qe'Ett 


0.2 


0.5 


0.6 


mass tuning 


0.2 


0.1 


0.2 


finite volume 


0.1 


0.3 


0.3 


relativistic correct. 


1.0 


1.0 


0.0 


operator matching 


4.1 


4.1 


0.1 


Total 


4.4 


4.6 


1.5 



6. making the coefficients ci depend on aAf;,; 

7. using (A.l) rather than (A. 7) for 5fq\ 

8. allowing for a 20% error in the scale f — Jtt (see 
appendix A for the relevant formulas). 

Fig. 9 summarizes results from these tests. We compare 
fs at the physical point with these modifications in place 
with results obtained with the basic ansatz. The latter 
corresponds to the left most data point in Fig. 9 and is the 
same as the magenta point in Fig. 8. The integers on the 
horizontal axis in Fig. 9 refer to the type of modification 
of the basic ansatz as enumerated above. One sees that 
the basic ansatz result is very stable. The decay constant 
fs changes by less than IMeV in all the tests undertaken. 

In Fig. 10 and Fig. 11 we show chiral/continuum extrap- 
olations of $s = /b^a/MbJ and ^s/^ both carried out 
and tested along similar lines as for Jb in Fig. 8 and Fig. 9. 
The x^/dof for the two extrapolations were 0.59 and 0.48 
respectively. 

The physical point results in Fig. 8 and Fig. 10 show sta- 
tistical, extrapolation and rp^^'' errors whereas in Fig. 11 
statistical and extrapolation errors are included. In the 
next section we will discuss additional systematic uncer- 
tainties inherent in our decay constant determinations, 
the dominant such error coming from higher order oper- 
ator matching. 

VI. RESULTS 

Table V gives the error budget for /s, Jb^ and /bJ Jb- 

3/2 

The first four entries, "statistical", "scale ", "contin- 
uum extrapolation" and "chiral extrapolation & gB'Bn" 
are all part of the errors emerging automatically from the 
chiral/continuum extapolations and shown at the physi- 
cal points in Fig. 8 and Fig. 10. Their individual contribu- 
tions were separated out using the methods of reference 
[isj (see eq.(30) and (31) of that article). The remaining 
four entries in Table V, "mass tuning" , "finite volume" , 



"relativistic corrections" and "operator matching" are 
additional systematic errors affecting our calculations. 
Sensitivity to the strange quark mass can be estimated 
by comparing results for valence quarks masses arus and 
anil . Similarly effects of mistuning of aMb can be investi- 
gated using older NRQCD/HISQ decay constant results 
(see [l^) covering a range of aMb values. Those calcu- 
lations were done before proper retuning of the 6-quark 
mass and provide information on how the decay constants 
depend on aMb. For the "finite volume" uncertainty we 
take the same percentages as determined for the D and 
Ds meson decay constants in [2^ using finite volume chi- 
ral perturbation theory. Our heavy-light currents have 
been matched to full QCD through order asAgcD/Mb 
and corrections come in at order {h^Qco/MbY ~ 0.01. 
There are order agAQco/Mb corrections to the NRQCD 
action that are not included in our simulations. However, 
as discussed in Q , their effect on decay constants can be 
bounded to be at most ^1%. 

The O(cks) corrections to eqs.(I7][5|) are not known. The 
j{t},sub nonleading, so the most important high-order 
correction is in the coefficient of Jq°^ . To account for cor- 
rections at this level and beyond, we modify our data by 
multiplying the right-hand side of eq.(l7]) by an overall 
factor of 1 -|- o^Pq where we approximate w 0.1. We 
use two different p'^s, one for all coarse-lattice data and 
the other for all fine-lattice data. To be conservative, we 
take each to be 0(0.4), which is more than twice as large 
as the one-loop pa and also comparable to the estimates 
used in @: that is, we set each p'q = ± 0.4. The er- 
rors from these factors are combined in quadrature with 
the simulation errors in the currents, taking care to pre- 
serve the correlations caused by the fact that all course- 
lattice data has the same pg , as does all fine-lattice data. 
We then repeat the chiral/continuum extrapolations de- 
scribed in the previous section, this time applied to the 
modified data with enhanced errors. We use the differ- 
ence between the total extrapolation error obtained with 
and without higher order matching errors added to the 
data to estimate the operator matching errors for /b and 
Jb, ■ These are given as the last entries in Table V. 

Finally, we note that sea charm quarks are omitted in 
our simulations. However we expect their contributions 
to be small enough that the final total errors in Table V 
are unaffected. 

Our final decay constant results including all the errors 
discussed above are, 

fB = 0.191(9)GeV, (33) 

/s, = 0.227(10)GeV, (34) 

and 

^ = 1.188(18). (35) 

JB 

These numbers are in good agreement with HPQCD's 
previous NRQCD 6-quark/AsqTad light quark results [53, 
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HPQCDthis work '12 
using fgj from HISQ b 
directly from NRQCD b 



HPQCD NRQCD b '09 



Fermilab/MILC ' 1 1 



ETMC ' 1 1 (Nf=2) 



-e- 



e- 



120 140 160 180 200 
B Meson Decay Constant in MeV 



220 



HPQCD this work '12 

HPQCD HISQ b ' 1 1 
HPQCD NRQCD b '09 

Fermilab/MILC ' 1 1 

ETMC ' 1 1 (Nf=2) 



-B- 



-e- 



160 180 200 220 240 260 
Meson Decay Constant in MeV 



FIG. 12: Comparisons of results for fs from this article 
with previous HPQCD work and with results from the 
Fermilab/MILC ^ and ETM ^ collaborations. 



FIG. 13: Comparisons of results for /b^ from this article 
with previous HPQCD work [1, and with results from the 
Fermilab/MILC [3] and ETM ^ collaborations. 



however with improved total errors. Comparison plots 
are shown in the next section. 

The errors in /b and /b, are overwhelmingly domi- 
nated by the matching uncertainties. Without them, the 
total errors would be reduced to 4.6% — >■ 2.1% and 4.4% 
— ^ 1.6% for Jb and Jb^ respectively. Clearly a huge ad- 
vantage can be gained if one could develop a formalism 
that did not require operator matching. One major mo- 
tivation for designing the HISQ action Q was to come 
up with a quark action that could be used not only for 
accurate light quark physics, but also to simulate heavy 
quarks. It has been e niploye d already very succesfully 
for charmed quarks [l^, [20l - [2^ and HPQCD has recently 
also started work with amq > arUc [13 ■ The HISQ ac- 
tion allows for a relativistic treatment of heavy quarks 
which means that one does not have to resort to effec- 
tive theories. One important consequence is that decay 
constants can be determined from absolutely normalized 
currents. There is no need for operator matching. Fur- 
thermore it has been demonstrated that due to its high 
level of improvement the HISQ action can be used for 
heavy quarks up to about amq < 0.8 without leading 
to large discretization effects. Recently a succesful appli- 
cation of heavy HISQ quarks to B physics was achieved 
through a very accurate determination of the Bs meson 

decay constant, namely f^^^^^ — 0.225(4)GeV with er- 
rors of only 1.8% [23]. There is very good agreement 

between f'j^/'^'^^ and the NRQCD 6-quark resuh eq.([34l) 
of this article. This indicates that the very different sys- 
tematic errors in the two calculations are under control 
and properly accounted for in our error estimates. 



The HISQ 6-quark calculation of fs^ required going 
to very fine lattices including the MILC superfine and 
ultrafine ensembles with lattice spacings ~ O.OGfm and 
~ 0.045fm respectively. Repeating those calculations for 
the B meson with its light valence quark would be quite 
expensive and it will take some time before such calcula- 
tions become available. In the mean time we can combine 
j.(HiSQ) the result eg. (1551) of this article to extract 
a new and accurate Jb- One finds. 



fB 



AHISQ) 



fB = O.I89(4)GeV, (36) 



NBQCD 



which is in excellent agreement with (1331) . only more ac- 
curate by better than a factor of 2. Eq. ([36| is the most 
important result of this article for phenomenology. It 
also demonstrates the advantages of working with both 
HISQ and NRQCD &-quarks in parallel. In the future 
we plan to apply this combined approach to B and Bg 
semileptonic decay studies as well. 



VII. SUMMARY 

We have carried out new determinations of /s, Jb^ and 
JbJJb using NRQCD 6-quarks and HISQ hght valence 
quarks and improve on our previous calculations with 
AsqTad light quarks. Figures 12, 13 and 14 c omp ares 
our new results with HPQCD's older work [a, |23| and 
also with results from the Fermilab/MILC [13] and the 
ETM 24] collaborations. One finds overall consistency 
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HPQCDthis work '12 



HPQCD '09 



Fermilab/MILC ' 1 1 



ETMC ' 1 1 (Nf=2) 



-e- 



■e- 



-e- 



0.9 



1.1 

Ratio fgj, / fg 



1.2 



1.3 



fg from Lattice QCD: 

HPQCDthis work, '12 



fit 



from global fits: 



(Lunghi and Soni '11) 

no [sin(2 beta)] input 
no [B -> tau,nu] input 



120 140 160 180 200 220 
B Meson Decay Constant in MeV 



FIG. 14: Comparisons of results for fB^/fs from this article 
with previous HPQCD work [5|| and with results from the 
Fermilab/MILC jT,] and ETM ^ collaborations. 



FIG. 15: Comparisons with fg^^^ from global fits as given 
in[|. 



Appendix A: Partially Quenched ChPT Chiral Logs 



between the different lattice groups. Our most accurate 
determination of fs, ea. p6p . comes from combining the 
new ratio fsj fs, (l35l) . with a precise determination of 
fs^ based on HISQ 6-quarks This gives the most 

precise Jb available today with errors of just 2.1%. As 
mentioned in the introduction, accurate values for Jb are 
needed to compare with from global fits in Unitar- 

ity Triangle analyses. In Fig. 15 we compare the new 
accurate Jb with two examples of fg^^"^ determined by 
Lunghi and Soni With current errors the two fg^^'^ 
values are consistent with each other and with /s from 
Lattice QCD. In the future, once errors are reduced con- 
siderably, these kind of comparisons could become more 
interesting. 
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In this appendix we summarize partially quenched 
ChPT (PQChPT) expressions for the chiral logarithm 
terms S/b and SJb^, taken from the literature. We fol- 
low closely the notation of 1^ which we also adopted in 
our D —i' K semileptonic paper JjQj . We use "u" and "s" 
for sea and "g" and "g^" for valence light and strange 
quarks respectively. Furthermore mab is the mass of the 
pseudoscalar meson with quark/antiquark content a and 
b and m2 = ^{ml^ + 2ml^). For x = q or q, PQChPT 
gives. 



5fE 



327r2/2 



-2Ii{mxu) - h{mxs) 



where 



and 



with 



Ii (m) = m^log 



m 
A2 



(Al) 



(A2) 



DRi'^'\m;I) = ^R^'-'^im;I), (A3) 



(m^ — m^) 
[m? — rw^) 



I{m) 



(A4) 
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is the B*Btt coupling which is only known approxi- 
mately from experimental data in the D system. In our 
chiral/continuum extrapolations we treat this "constant" 
as one of the fit parameters and set its priors to a central 
value of 0.25 with width 0.10 (40%). The scale A is set to 
47r/, with / given by the physical pion decay constant. 
In the full QCD limit the partially quenched formulas 
simplify to, 



with, 



3.9^ 



327r2/2 



[-2J{mxu, A + Sxu) - J{rnxs, A + Sxs) 



(A7) 



SfE 



SfB 



327r2/2 



1 + 3.g2 
327r2/2 



-2/i(mA') - -/i(m^) 



(A5) 



3 1 

--Ii{m.^) - -/i(m,,) - IiirriK) 



(A6) 



Following we have also considered PQChPT logs that 
include hyperfine and flavor splittng effects. A modifica- 
tion of the terms proportional to 352 in (|Aip is required 
leading to. 



327r2/2 



[-2Ii{mxu) - h{mxs) 



J(m,A) = {m^-2A^) log 



F{l/x) 



A2 



tan 



,2A2-4A2F(m/A), 
(A8) 



\x\ < 1 



^^^^ log{x + Va;2 - 1) 



\x\>l 
(A9) 

A is the B* — B^ hyperfine splitting and 5xu and 5xs 
adjust for the fact that in some one-loop diagrams the 
internal S*/^ does not have the same fiavor as the exter- 
nal Bx- We have carried out chiral/continuum extrap- 
olations with both (A.l) and (A. 7). Differences in the 
final values at the physical point serve as a measure of 
systematic errors coming from our extrapolation ansatz. 
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